Photocurrent measurements of supercollision cooling in graphene 
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The cooling of hot electrons in graphene is the critical process underlying the operation of exciting 
new graphene-based optoelectronic and plasmonic devices, but the nature of this cooling is contro- 
versial. We extract the hot electron cooling rate near the Fermi level by using graphene as novel 
photothermal thermometer that measures the electron temperature (T(t)) as it cools dynamically. 
We find the photocurrent generated from graphene p — n junctions is well described by the energy 
dissipation rate CdT/dt — — A(T 3 — Xf ), where the heat capacity is C — olT and Ti is the base 
lattice temperature. These results are in disagreement with predictions of electron-phonon emission 
in a disorder-free graphene system, but in excellent quantitative agreement with recent predictions 
of a disorder-enhanced supercollision (SC) cooling mechanism. We find that the SC model provides 
a complete and unified picture of energy loss near the Fermi level over the wide range of electronic 
(15 to ~3000 K) and lattice (10 to 295 K) temperatures investigated. 



How does an excited electron lose its energy? This is a 
central problem in fields ranging from condensed matter 
to particle physics. One key pathway is the emission of 
massless bosons such as photons or phonons. However, 
momentum must be conserved and the phase space avail- 
able for such emissions can be dramatically restricted. 
For example, an electron moving through free space can- 
not emit a photon without transferring momentum to a 
third body. 
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FIG. 1. Hot electron cooling by acoustic phonons. (a) 

Momentum conservation restricts cooling of hot carriers near 
the Fermi level (green) to low energy (< 4 meV) acoustic 
phonon emission (black arrows, scale exaggerated), (b) In 
a supercollision (SC) transition the momentum restrictions is 
relaxed by the lattice disorder (q d ), enabling faster cooling by 
emission of high energy (~ ksT) acoustic phonons. Solving 
the SC rate law dT/dt = -(A/a)(T 3 - T z 3 )/T, we plot the 
predicted cooling of the graphene hot electron temperature 
T(t,Ti) — Ti (log-scale). The thermal decay changes from 
inverse to exponential with increasing lattice temperature. 



The cooling of hot electrons in graphene presents an- 
other interesting case. Here, hot electrons (e - ) move 
at a constant speed on a conical energy-momentum 
surface, and dissipate heat by phonon emission. The 
optic phonon energies in graphene are unusually high, 
huj op > 200 meV and mediate cooling for only very 
hot electrons. pQ For electrons with energy below hw opi 
acoustic phonon emission is the dominant cooling path- 
way. However, these phonons move with the much slower 
sound velocity v s <^vp- [2,3 As shown in Fig la, this 
velocity mismatch, combined with momentum conserva- 
tion, greatly restricts the energy (E ac ) of emitted acous- 
tic phonons to AE ac /kBT < 2v s /vp ~ 0.04 , where ksT 
is the typical energy of a hot electron. More than ~40 
acoustic phonons to cool a hot electron to just half an 
initial energy of 0.2 eV.[2j[4] This inefficient process cre- 
ates a cooling bottleneck, with calculated cooling times 
exceeding 300 ps.[3j[5] 

Alternatively, a recent theory by Song et al. predicts 
that disorder effectively relaxes the momentum conser- 
vation constraint, enabling the emission of large energy 
(ksT) and momentum (kBT/hv s ) acoustic phonons, as 
shown in Fig. lb. This mechanism is called supercollision 
(SC) cooling, and the theory predicts relaxation times of 
1-10 ps, orders of magnitude faster than the disorder- free 
model. [6] 

Here, we perform the first experiments to directly test 
the conflicting predictions of hot electron models. [TJ [3j 
HJ [6l [7] We use the photothermal effect in graphene as a 
novel quantitative probe of hot carrier cooling near the 
Fermi level. We find excellent agreement with the predic- 
tions of the supercollision model, showing that disorder 
effectively relaxes momentum conservation and leads to 
very rapid electron cooling. Using the cooling rates ex- 
tracted, we directly determine the hot electron temper- 
ature in graphene, which is of central importance both 
to graphene's fundamental physics and for its use in a 
variety of electronic and optoelectronic applications such 
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FIG. 2. Photocurret setup, a time-resolved graphene thermometer, (a) The temperature T(t) of thermalized hot 
electrons and holes in a graphene p — n junction cool at a characteristic rate r _1 . The elevated junction temperature drives 
the collected thermoelectric current, i(t). (b) Experimental setup where we collect PC from graphene as a function of TJ, laser 
power (P) and two-pulse time delay (td)- Graphene (green, false color) is excited using a focused CW or pulsed laser. The 
device SEM shows an overlay of a spatial PC map with peaks for the graphene p — n (0.7 nA, red) and n — p junctions (-0.6 nA, 
blue), (c) Pulsed excitation PC map, plotting electrons collected (Qi/e) vs. applied gate voltages. Tuning the electrostatic 
gates (Vtg, and Vbg) show six PC regions (dotted lines). 

as photodetectors and bolometers. [5HTT] 

The energy relaxation rate of a hot electron gas, ^ = 
C ^ = Pi n — H is determined by the heat loss rate H and 
the heat capacity C,[3] where Pi n is the incident power 
delivered to the electrons. For a degenerate electron gas 
with heat capacity C = aT and T > Tiv s ^ (typically 
only 5-10 K), the SC mechanism shown in Fig. lb pre- 
dicts Hsc — A(T 3 ; — XJ 3 ), where A is rate coefficient, k F is 
the Fermi momentum and T\ is the lattice temperature. [6 
For comparison, the conventional momentum conserving 
model (Fig. la) gives H = A'T A (T - T t ) for E F <C k B T 
or H = A"(T - Ti) for E F > k B T, where E F is the 
Fermi energy. [T]-[3] 

Under steady-state conditions where H = P ini the SC 
model predicts the following temperature scaling with 
input power: 

T = (Pin/A) 1 / 3 , T > Ti 

T = 2] + 3^fp T-T i «T i . (1) 

If instead we deliver a short impulse of energy P in = 
F in 5(t) to the system, the electron gas is heated to an 
initial temperature of T Q = \/Tf + 2F in /a, where F in 
is the remaining deposited energy after the initial opti- 
cal phonon heat dissipation. [12] The subsequent decay of 



transient electron gas temperature T(t) is governed by 
dT/dt — —Hsc/C with solutions: 

1 + t/T 

T(t) = T, + (T - T,)e^ , T(t) - T, < T, (2) 

where r" 1 = (A/a)T and rf 1 = (3A/a)T t are charac- 
teristic hot electron cooling rates. The full solution for 
T(t,Ti) is plotted in Fig. 1 using a rate coefficient A/a 
we later determine as 5.5 x 10 8 K -1 s -1 . With increasing 
7), the thermal cooling, T(t, T{)1 changes from inverse to 
exponential in time. 

To experimentally test the above predictions, we lo- 
cally heat a graphene p — n junction with laser and use 
photocurrent generated as a thermometer of either the 
steady-state (Tew) or transient (T(t)) hot electron tem- 
peratures. When heating graphene using 180 fs-long light 
pulses, we assume that only a fraction 7 = F in /F of the 
total incident laser pulse energy (F) is retained in the 
hot electron gas created. This thermalized distribution 
is characterized by an initial temperature T , and cools 
dynamically at a rate r _1 (see Fig. 2a). Similarly un- 
der continuous wave (CW) illumination, only a fraction 
7 = Pin/P of the total incident laser power P is coupled 
into the electron gas, maintaining a steady-state temper- 
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at lire. 

In Fig. 2b, we show a schematic of the photocurrent 
measurement setup and single-layer graphene p — n junc- 
tion photodetector device. The junctions are created by 
globally p (or n) doping the graphene sheet with an elec- 
trostatic back gate (BG), and locally n (or p) doping 
through a top gate (TG) . [9] H3] We further overlay a spa- 
tial PC map on our device SEM image: Positive (red) or 
negative (blue) PC peaks are measured as we raster scan 
a 1.5 /am diameter laser spot over the p — n and n — p 
junctions, labeled. Data is collected at T/=10 K unless 
otherwise indicated. 

Figure 2c plots the charge (Qi) collected per excita- 
tion pulse using a laser repetition rate (/) of 76.1 MHz. 
Each pulse induces a time-dependent photocurrent re- 
sponse i(t). We measure the resulting integrated charge 
Qi — J i(t)dt, or the average PC given by Qif. As 
the applied gate potentials are tuned, a sixfold pattern 
of alternating-sign photocurrent emerges, corresponding 
to p — n, p — p + , p + — p, n — p, n — n + , and n + — n 
junctions. A similar pattern in also observed under CW 
excitation (see supplementary materials). It was recently 
shown in graphene that such six-fold PC patterns indi- 
cate electron-hole separation occurring by a thermoelec- 
tric process (illustrated in Fig. 2a). [TU US]). We can 
therefore use the measured thermoelectric current given 
by US HE1: 



i(t) = PT(t){T{t)-T l ). 



(3) 



to extract the hot electron temperature for both CW and 
pulsed excitation conditions. Here /3 is proportional to 
the Seebeck coefficient and is theoretically predicted to 
be ~ 2 pA/K 2 (see supplementary information). 

In Fig. 3a, we plot the PC collected at a p — n junction 
(Vtg=2 V, Vbg=-15 V) as a function of CW laser power. 
The photocurrent is sublinear and is accurately fit by a 
power law, growing as lew ~ P°- 65±0 02 . Figure 3b show 
an identical measurement using pulsed excitation. Com- 
paring the two excitation techniques at identical laser 
powers, the amplitude of PC generated under pulsed ex- 
citation is at least 10 times smaller than the CW case, 
as was reported previously. [17] Similar to the CW pho- 
tocurrent, we find the current grows with a power law, 
but this time with: I lp ex p°- 50 ± - 03 (see Fig. 3b inset). 

To compare the above power laws extracted against 
the SC model we combine equation (1) with the thermo- 
electric model to predict a CW PC power dependence 
of: 



lew = PT* W = (3(P zn /A) 2 ^ 



lew 



3ATi 



T>T Z 

T;«T;. (4) 



The fitted power laws are I CW oc pO-65±o.02 for ^ = 1Q 
K and I C w oc P for T% = 295 K (Fig. 3c), in excellent 
accord with the SC model. Loss rates associated with 
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FIG. 3. PC response obeys SC power laws, (a) PC 

generated under CW excitation at a graphene p — n junc- 
tion scales as p°- 65±0 - 02 (black line), (inset) Correspond- 
ing electron temperatures scale linearly with \fP. (b) PC 
(Qif) and corresponding initial temperatures (T Q ) vs. pulsed 
laser power. Black line, power law fit of p - 50 ^ - 03 . (inset) 
Same plot converted to electrons collected per pulse (Qi/e) 
vs. photon fluence (square root scale). (c) CW power depen- 
dence from n — p junction becomes increasingly linear over 
the Ti = 10-295 K range. We estimate the cross-over current 
/ c , and use it to determine (5. 
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FIG. 4. Extracting the hot electron relaxation 
time, (a) Collected two-pulse photocurrent (Q 12 (td)f) re- 
sponse at selected incident photon fluences (in units of xlO 14 
photons/cm 2 ) (b) As shown, the decay of the TPC signal 
(AQ12 (td)/, normalized) is closely linearized when plotted 
on an inverse scale, (c) TPC peak amplitude (square root 
fit, dashed line) and r (inverse root fit, black line) vs. laser 
fluence, F. (inset) The product of the data points yields a 
constant AQi2(0)/r o = 160 =b 13 ps-pA (dashed red line) or 
an electron cooling rate of A/a of ~ 5.5 x 10 8 K~ 1 s~ 1 . 



other proposed momentum-conserving models H ex T or 
H oc T 5 predict powers that are well outside of the error 
bars of the measured exponent. [3] 

The pulsed excitation power dependence can also be 
predicted using the SC model temperature T(t) from 
equation (2): 



-F 

Jo 



i(t)dt = (3(a/A)T Q . 



(5) 



The total current Qif collected is thus linearly related 
to the initial hot electron temperature, which from above 
is T Q = ^2jF/a. Hence the resulting PC should scale 
as ft oc y/F, in excellent agreement with the data fits 
shown in Fig. 3b. 

The coefficient j3 can also be extracted from Fig. 3 by 
finding the CW cross-over current (J c ) where the lew 
power dependence transitions from ~ Ps to linear in 
P. We find f3 = ^I c /Tf (see supplementary materials). 
Figure 3c shows I c occurs at higher powers as the base 
temperature is warmed to 295 K. We read-off cross-over 
currents of -85 pA (T t =10 K), 1.2 nA (35 K) and 2.4 
nA (55 K) and calculate a mean /3 of 1.1 pA/K 2 for 
our device. In Fig. 3a (inset) we use this f3 to plot 
the graphene p — n junction temperature vs. incident 
power. These direct measurements of the graphene 
electron temperature are important for the design and 
feasibility of graphene device exploiting electron thermal 
gradients. [Hill] 

We have shown the SC model coupled with the 
thermoelectric effect predicts the functional form of the 
CW and pulsed PC measurements. However, these PC 
power trends do not directly measure the timescales 
for electron cooling, nor the associated hot electron 
cooling rate A/ a needed to quantitatively compare to 
the SC model and determine the absolute hot electron 
temperature in graphene. 

We use a time-dependent two-pulse excitation tech- 
nique to measure T(t) in graphene and extract the cool- 
ing rate. The experimental setup is outlined schemati- 
cally in Fig. 2b. The first pulse creates high energy e — h 
pairs at the graphene p — n junction, which rapidly ther- 
malize and cool to a temperature T Q on a rapid < 300 fs 
timescale associated with optic phonon emission. [I2j [18] 
The resulting distribution of hot electrons cools from T Q 
to a transient temperature T(td) by acoustic phonons at 
a characteristic rate r" 1 . At the pulse delay time a 
second collinear pulse of equal intensity is absorbed, heat- 
ing the electron gas to y/T Q + T(tdj. The resulting total 
charge Q 12 (td) collected will then vary with time-delay 
as the transient p — n temperature (T(td)) cools. 

In Fig. 4a the collected transient photocurrent (TPC) 
signal, Qi2(td)f is plotted for selected photon fluences F. 
As td — >• ps, the magnitude of PC collected is greatly 
diminished because of the sublinear dependence of the 
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PC on laser power (Fig. 3b). Analogous time-dependent 
reductions in PC have recently been recently reported 
for graphene based devices. [TJ [19] In Fig. 4b, we plot 
AQi2(td) on a normalized reciprocal scale. The TPC 
decay kinetics are not exponential, but instead show a 
striking resemblance to the 1/td thermal decay predicted 
by the SC model in equation (2). 

To quantitatively interpret these results we integrate 
the time-dependent photothermal effect from equation 
(5) piecewise about t dl giving: 

J t i(t — td, s/Tq + T(td) 2 )dt, where T Q is the initial tem- 
perature created by each pulse independently. Solving 
using the thermal decay in equation (2) we obtain: 

AQ 12 (t d ) = (3(a/A) (t o + T(t d ) - ^ + T(0) . 

(6) 

This resulting TPC response function for Qi 2 {td) is pro- 
portional to the transient temperature T{td). Fig. 4b 
(red lines) shows this analytic solution fits our data 
well, requiring only two parameters; the amplitude Qi — 
f3{a/A)T Q and the thermal decay rate, r" 1 = AT Q /a. 
This functional form further fits our data well for a wide 
range excitation wavelengths investigated (0.8 to 1.55 eV, 
data not shown). 

The TPC decay in Figures 4a and b becomes markedly 
faster with increasing photon fluence F. This strong flu- 
ence dependence is captured by the SC model thermal 
decay in equation (2), which states r Q = a/AT" 1 , or 
equivalently that r Q scales with 1/y/F. Plotting the ex- 
tracted fit parameters in Fig. 4c, we show hot electron 
cooling time (r , orange circles) decrease from 6.3 to 1.3 
ps, closely scaling with 1/y/F (solid line fit), as predicted. 
The transient amplitude AQi 2 {td = 0) also scales non- 
linearly as y/F (dotted line in Fig. 4d) up to a maximum 
F of 3 x 10 14 photons-cm -2 where the transient response 
saturates. Their product AQi 2 fr shown in Fig. 4d {in- 
set) is approximately constant. 

In the SC model, Q\t q = f3{a/A) 2 measures the fun- 
damental cooling rate coefficient A/ a. Using /3=1.1 
pA/K 2 found earlier, we find a SC cooling rate of A/a = 
(5.5 ±0.4) x 10 8 K -1 s -1 . 

Theoretical estimates of the SC cooling are given in 
Song et al. asl6 : A = ^ A¥ = I A¥ where the 

& L J a 7r z kpl n 3 kpl n 

electron-phonon coupling strength is A = ~2 ^^fy [£] 
Using estimates for the deformation potential, D = 10 — 
30 eV, Ep=0.1 eV and a mean free path of kpl = 10, 
this theory predicts: A/ a = 1 • 10 8 - 1 x 10 9 K" 1 s" 1 . 
(The range comes from the uncertainly in D). The best 
match to our experiments indicate D — 12 — 18 eV, well 
within the expected range. 

With our cooling rate coefficient now extracted, we 
now plot in Fig. 3b {inset), the initial temperature, T Q 
for the thermalized electron gas. T Q can exceed 1000 
K, an order magnitude higher than in the CW case. 
Once heated, our data predicts hot electrons cool with 
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FIG. 5. SC model predicts TPC dependence on T z . 

At constant incident power, the TPC response varies consid- 
erably upon warming to 295 K. Using the SC temperature 
response in Fig. 1, we calculate the predicted TPC response 
AQi2(£d) with no free parameters (gray lines). 



a relaxation time varying inversely with T , as r Q = 
{{A/a)^)- 1 =1.8ns/T [K]. 

Since all the parameters in the model have been de- 
termined, the SC model predicts the lattice temperature 
dependence of the transient electron temperature with 
no free parameters. Figure 5 plots TPC data for differ- 
ent base lattice temperatures for a constant laser pho- 
ton fluence of 1.1 x 10 14 photons/cm 2 , corresponding to 
T Q = 1250 K. Upon warming the lattice to room temper- 
ature, the amplitude of TPC signal shrinks by a factor of 
~3, and the kinetics exhibits a dramatic shift toward 
a rapidly decaying exponential function. To compare 
with theory, we use the analytic SC model solutions for 
T(t, T{) — T\ plotted in Fig. 1 to numerically solve for the 
TPC response, AQi2(td, T/). With no adjustable param- 
eters, the SC model curves in Fig. 5 accurately predicts 
both (i) the amplitudes, and (ii) strongly varying func- 
tional decay observed. Recent graphene time-resolved 
THz experiments report a similar change in decay kinet- 
ics with increasing TJ, as observed in Fig. 5 for our TPC 
data. [SHED 

The above results definitively show that the SC model 
gives an excellent quantitative description of both the 
CW and pulsed PC experiments. As a last demonstration 
of this connection, we connect the disparate magnitudes 
of the PC measured in the CW (Fig. 3a) and pulsed (Fig. 
3b) excitation. The ratio, lew /hp is predicted by the 
SC model using straightforward algebra to be 1/ \J^jr or 
equivalently ~ \.2y/T~ (see supplementary information). 
For the data shown in Fig. 3b, T Q ranging from 250 K 
to 3500 K, giving corresponding lew /hp ratio ranging 
from 8 to 18. This range is in excellent accord with the 
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10 to 20 range observed in Fig. 3, and provides an in- 
dependent check that CW and pulsed PC experiments 
can be explained by the same fundamental underlying 
physics of SC hot electron cooling. 

In summary, we have introduced a quantitative frame- 
work for interpreting CW, one and two-pulse PC exper- 
iments as measurements of hot electron cooling of elec- 
trons near the Fermi energy. Over a broad range of elec- 
tron (20-3000 K) and lattice (10-295 K) temperatures, we 
find the electron gas heat loss rate is Hsc = A(T 3 — XJ 3 ) 
with a rate coefficient A/a = 5.5 x 10 8 K _1 s _1 for our 
device. At low lattice temperature the associated cooling 
time given by r Q = [(A/a)T ] _1 = 1.8 ns/T [K]. These 
cooling times are much faster than those predicted by 
acoustic phonon emission but are in excellent agreement 
with disorder-assisted supercollision cooling. The cooling 
rates extracted directly determine the graphene electron 
temperature, which is of central importance in design- 
ing graphene terahertz plasmonic devices, photodetectors 
and bolometers. 



METHODS AND MATERIALS: 

Single-layer graphene on copper foil is grown using the 
Chemical Vapor Deposition (CVD) method [22]. Micro- 
Raman was used to confirm the growth of large-grain 
single-layer graphene with no visible D peak, indica- 
tive of high quality growth. Graphene was transferred 
using the lift-off technique onto a 300 nm Si02 layer 
grown on top of a silicon wafer which serves as the global 
back-gate (BG). The large-grain growth graphene is di- 
vided into 30x50 jam stripes using photolithography fol- 
lowed by oxygen plasma etching. Electrode pads of ti- 
tanium/gold (3nm/150nm) are deposited along graphene 
stripes with variable source-drain distances of 10 or 20 
/im. A good dielectric separation with the top gate is 
achieved with 10 nm of Si02 by electron beam depo- 
sition, followed by Hf02 atomic layer deposition. Fi- 
nally, an optically translucent top gate of titanium/gold 
(2nm/20nm) is deposited along the center of the source- 
drain gap with a width of 6 jam. The device is mounted 
in an Oxford HI- RES liquid helium cryostat. The CVD 
graphene photodetector device had a characteristic high 
mobility of ~8,000 cm 2 V _1 s _1 with centrally located 
Dirac points in conductance sweeps (see supplementary 
information). 

Light was generated by a Coherent MIRA oscillator 
that was externally compressed with via a prism-pair 
line. Autocorrelation measurements at the cryostat po- 
sition yield beams centered at 1.25 eV and show a 180 
fs FWHM pulse duration. For TPC measurements the 
beam paths were cross-polarized to suppress pulse inter- 
ference effects. After a mechanical delay stage, the two 
beams are aligned in a collinear geometry at a beamsplit- 
ter (BS) and scanning mirror. They are coupled into the 



microscope (Olympus BX-51) through a 50XIR Olympus 
objective with cover glass correction and piezo scanning 
mirror (SM, Fig. 2b). The TPC response is collected 
as function of time-delay (td) at 1 kHz beam modulation 
via lock-in and current amplifiers. 
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